Many-body theory of the quantum mirage 



o 
o 

(N 



i 

CO 

c3 



I 

C 

O 

o 



> 
(N 
oo 
o 



C3 



T3 

c 

O 
o 



X 
S3 



A. A. Aligia 

Comision Nacional de Energia Atomica, 
Centro Atomico Bariloche and Instituto Balseiro, 
8400 S.C. de Bariloche, Argentina. 
(Received February 1, 2008) 

In recent scanning tunneling microscopy experiments, confinement in an elliptical corral has been 
used to project the Kondo effect from one focus to the other one. I solve the Anderson model at 
arbitrary temperatures, for an impurity hybridized with eigenstates of an elliptical corral, each of 
which has a resonant level width 5. This width is crucial. If S < 20 meV, the Kondo peak disappears, 
while if S > 80 meV, the mirage disappears. For particular conditions, a stronger mirage with the 
impurity out of the foci is predicted. 



The recent advances in nanofabrication has led to fas- 
cinating experiments. In one of them, the electrons on 
the Cu(lll) surface are confined in a quantum corral as- 
sembled by depositing other atoms on the boundary of 
an ellipse. When a Co atom is positioned at one focus 
of the ellipse, the corresponding Kondo signature is ob- 
served by scanning tunneling microscopy (STM) not only 
at the atom, but also at the other focus. jjJD Another 
remarkable feature of the experiments is that eigenstates 
of a two-dimensional free electron gas confined into the 
ellipse are clearly displayed. The average separation be- 
tween the corresponding eigenvalues is ~ 10 meV, and 
one might expect that the confinement is good enough 
such that the level width of each eigenstate 6 < 10 meV. 
On the other hand, a Kondo peak in the impurity spectral 
density of states is in general observed when the impurity 
hybridizes with a continuum of conduction states. Can 
the Kondo feature exist if these states are discrete? From 
recent work in mesoscopic systems, we know that the an- 
swer is positive if the Kondo temperature Tk <^ d, where 
d is the average separation of the levels which have a sig- 
nificant hybridization with the impurity. [0] However, in 
the experiments Tk ~ 5 meV, while d ~ 100 meV. 

In this Letter, I calculate the impurity and conduction 
density of states as a function of 5. This requires an 
explicit calculation of the many-body effects, which has 
not been done in previous theories of the mirage effect. 
P,|]-0 I find that only for values of 5 large enough to lead 
to a well defined Kondo peak, but low enough to allow 
coherence between both foci, the theory is consistent with 
experiment. From the insight gained on the decoherence 
effects, I propose a slight modification of the experiment, 
which should result in an enhanced mirage effect, with 
the impurity out of the foci. 

To describe the essential features of the experiment, I 
use the following Anderson model: 

H = ^2 e i c ] a c j° + E d^2 d v d ° + ud \dfd\di 

ja a 

+XVi( c 5«A+H.c.)+ff'. (i) 



Here c- and d\ create an electron on the j th eigenstate 
of a hard wall elliptic corral and the impurity respec- 
tively. Roughness of the boundaries of the corral on a 
length scale smaller than the Fermi wave length 2tt /hp ~ 
30A does not affect much the wave functions. Jl],§ H' de- 
scribes the hopping of each corral state with a continuum 
of conduction states outside the corral, with the same 
symmetry. If the hopping and spectral density of these 
states is constant, and if the effective interaction between 
corral states can be neglected, the effect of H' is simply 
to introduce a width 8 in the unperturbed (Vj — 0) Mat- 
subara Green function of each state: 



G°Jiuj n ) = ((c jcr ;c] (J ))i LJn = [iuj n - e + i8sgn(u) n )]~ 



(2) 



I take this form with S independent of j for simplicity. Vj 
should be proportional to the (real) j th normalized wave 
function tpj(r) of the corral at the impurity position Ri. 
I take: 



Vj = 25 meV Vabifj(Ri) max I 1 



eV 







(3) 



a (b) is the semimajor (semiminor) axis of the ellipse, 
and £p is the Fermi energy. The last factor produces a 
slow decrease of Vj with energy and a smooth cutoff 1 
eV above Ep. It leads to a more symmetrical line shape, 
equilibrating approximately the weight of the Vj at both 
sides of ef and does not affect the physics. The energy 
prefactor was chosen to lead to the experimental width 
of the Kondo feature for S = 50 meV. If the Cu surface is 
represented by a tight binding model with one orbital per 
site, M this amounts to a hopping t' ~ 0.67 eV between 
the impurity and the atom below it (or t' /n if the impu- 
rity has n nearest neighbors on the surface). || This value 
is reasonable for a hopping between 4s and 3d electrons. 
The rapid decay of this hopping with distance allows to 
neglect hopping of the impurity with other atoms. Since 
in addition the density of s and p states at Ef is two times 
larger in the clean surface than in the bulk, |^| one can 
safely neglect the coupling of the latter with the impurity. 
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The eigenstates of free electrons inside an ellipse 
were obtained solving the matrix Hamiltonian in the 
basis of eigenstates for a circle, deformed to fit into 
the ellipse. p0| These basis states are proportional to 
J k { lknP ){e* k1r ± e- lfee ), where p = [(x/a) 2 + (y/6) 2 ] 1 / 2 , 
8 = arctan[ya/(a!&)], J k is the k th Bessel function and 
7fe„ is its n th zero. We retained the 100 lowest states 
in each of the four symmetry sectors of the group C% v . 
From the solution with an enlarged basis with a total 
of 1600 states, the error in the Ej near the Fermi en- 
ergy is estimated ~ 10 -6 eV. Recently, this problem was 
solved analytically. || We take an ellipse with eccentric- 
ity e = 1/2, a = 71.3 A (b = y/3a/2), and the effective 
mass was adjusted so that £42 = Sf- 

The particular structure of the non-interacting prob- 
lem introduces technical difficulties in the many-body 
problem. In addition, the fact that the temperature of 
the experiment T=4K<C Tk — 53 K, seems to invali- 
date the use of the non-crossing approximation, since it 
works for T > Tk and violates Fermi liquid relations 
at low T. Q,|ll[| I obtain the retarded Green function 
for the impurity Gd(ui) — ((d a ; d\))u using fmite-T per- 
turbation theory up to second order in U and analytic 
continuation: || G^ 1 = (G°) _1 - S, with G° d (iuj n ) = 
[iw„ -E d - Ej V?GQ(iu n )]-\ E d = E d + U(4d a ), and: 

E(tW,) = -{UT f Gd(i"l ~ lVrn)G d (iLU n )G d (lLU n + iv m ), 

(4) 

with u>„ = (2n + l)irT and v m = 2nnrT. For cases stud- 
ied before (with a simple band structure), comparison 
with calculations using Wilson's renormalization group 
H shows that the approximation is qualitatively cor- 
rect for all values of U. The method has proven to be 
quantitatively correct even for U = 2.57rA, where A is 
the resonant level width. [|l2| In the present case, the 
particular features of the non-interacting system renders 
an estimate of A difficult, but in absence of the corral 
A ~ 0.2 eV. @ I have chosen U = 1 eV. This value 
is enough to lead to a strong Kondo peak, remaining 
within the range of validity of the approximation. Ed 
was adjusted in order that the Hartree-Fock effective d 
level Ed falls near Ep, as suggested by first-principles 
calculations, |7],[l4| with a small shift of -22 meV to dis- 
place the Kondo feature to the experimental position for 
S = 50 meV. The resulting impurity spectral function 
Pd{ijj) = — Ini{Gd(w)}/7r is represented in Fig. 1 (a) 
for S — 50 meV. The shifts in the broad peaks from 
E d ~ -0.5 eV and E d + U ~ 0.5 eV are due to the 
particular structure of Ej and Vj . The only peak of inter- 
est here is the central one. Its temperature dependence 
is shown in Fig. 1 (b). The peak looses intensity and 
broadens as the temperature is increased. This is a clear 
signature of its many-body nature. For T = Tk, the 



maximum of the peak with respect to the background is 
reduced to roughly half its value at T = 0. 

Once Gd(u)) is known, the local conduction electron re- 
tarded Green function G c (r, u>) = ((i/)(r); ip(r)^))ui with 
ip(r) — Ej¥'i( r ) c j: can be obtained using the follow- 
ing exact relationship (simplified for real <fj(r)) obtained 
from the equations of motion: 

G c {r,u>) =^^ 2 (r)G°H + G d H{5(r,c)} 2 , (5) 
3 

with 

S(r,u)=J2<Pi(r)V j G° j (co). (6) 
3 
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FIG. 1. (a) Impurity spectral density as a function of en- 
ergy, (b) Temperature dependence of the central peak. 

The differential tunneling conductance dl/dV is essen- 
tially proportional to the local conduction spectral den- 
sity p c (r,u) = — Im{ G c (r, u/)}/ir. At or very near B4 
there might be a small deviation due to direct tunnel- 
ing between the STM tip and the impurity, which affects 
the line shape. Jll],[l3| This direct contribution cannot be 
transmitted to the mirage and for the sake of clarity, I 
neglect it here. In addition, the similarity between the 
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observed line shape at both foci suggests that this effect 
can be neglected. 

For the case of ipj(r)Vj independent of j and constant 
density of £j (as in the case of realistic models for 4f 
impurities in a metallic bulk |3j with r at the impurity 
site), the sum S(r, w) can be evaluated easily, and in usual 
cases in which the distribution of Sj is much broader than 
the peak in the resulting S lies near the imaginary 

axis for w near the peak, and S 2 can be approximated 
by a negative real constant. Then, the change in p c (r, u>) 
after addition of the impurity Ap c = — lm{GdS 2 }/ir, is 
proportional to -pd{u) near its peak. Physically, the level 
repulsion with the impurity states causes a depression of 
p c (r,u>). However, the above mentioned hypothesis are 
not valid in the present case, and although this picture 
remains qualitatively correct for large 5, the situation is 
different for small 8. 



The fact that S is dominated by the term with j = 42, 
means that essentially only this state "feels" the presence 
of the impurity at Ri and leads to a large mirage effect, 
since this state transmits coherently the information to 
the other focus (see Fig. 3 (a)). However, in the exper- 
iment the intensity at the empty focus is nearly 8 times 
smaller. As 8 increases, other conduction states increase 
their contribution to S and p c . By symmetry, all of them 
are even under reflection through the major axis, but 
they can have any parity under reflection through the 
minor axis a. While 9542(7") is even under a, the other 
states nearest to Sf with important hybridization with 
the impurity (\Vj \ > | V42 1/2), namely 32, 35 and 51, are 
all odd. These states contribute with the same sign as 42 
to the imaginary part of S at the impurity focus Ri (see 
Eqs. (||) and @), but with opposite sign at the empty 
focus aRi = —Ri- In other words, at the position of the 
mirage there is a negative interference of the states with 
opposite parity which tends to destroy coherence. For 
8 = 0.1 eV, the mirage is lost. 




FIG. 2. Impurity spectral density as a function of energy 
for several values of 8. 

In Fig. 2 we show the dependence of Pd{^) with 8 at 
T = 0. If 8 is very small, there are two narrow peaks, 
which correspond to two discrete energy levels for 6 — > 0. 
This agrees with the results of Thimm et al. for the case 
when one of the discrete states coincides with ep (the 
figure remains practically unchanged if Ef is increased to 
£43). |Q In addition, for 8 = 1 mcV, p c (r,uj) increases at 
the position of these two peaks after adding the impurity, 
and the shape of Ap c (u>), although smoother than that 
of p c (tjj) is also in strong disagreement with experiment. 
The structure with two peaks in Pd(w) disappears for 
5 ~ 18 meV, but a clear sharp peak does not develop 
until 8 ~ 30 meV. In spite of this, as shown in Fig. 3 
(a), Ap c (Ri,ui) for 8 = 20 meV is roughly consistent 
with experiment. This is easily understood noting that 
for small 6, S is dominated by the term proportional to 
G42 — (w — £42)^ , and lm{GdS 2 } becomes narrower 
than Im Gd ■ 
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FIG. 3. Change in the local conduction density of states 
at the impurity site (full line) and at the other focus (dashed 
line) for two values of S. 

To confirm this physical picture, I have repeated the 
calculation for 8 = 0.1 eV, setting Vj = for j =32, 35 
and 51. The Kondo feature becomes significantly nar- 
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rower and the mirage is restored with an intensity ratio 
1/2. From results similar to those shown in Figs. 3 and 
4 as a function of 5, the experimental situation seems to 
correspond to 8 ~ 40 meV. If the direct tunneling be- 
tween tip and impurity is indeed important, 8 might be 
smaller. I have verified that S approaches an imaginary 
constant for 8 > 50 meV, and that (as a consequence) 
Ap c (w) is quite similar to Pd(u) inverted for 8 = 50 meV 
and all temperatures displayed in Fig. 1 (b). 

In Fig. 4, I show the spatial dependence of p c and 
Ap c . The former displays some features of the state 43 
and is similar to the experimental topograph (Fig. 2 (c) 
of Ref. 0). However, it is more sensitive to ui and ep 
than A,o c , and to compare with experiment it is necessary 
to perform an integral over u>. Q In the difference plot 
Ap c (r), the contribution of states which do not hybridize 
with the impurity (like 43) is removed, and as in the 
experiment, Ap c (r) displays essentially minus the density 
of the state 42. One effect of increasing 8 is that the 
minimum of Ap c (r), which (as the maximum of |^42(?")| 2 ) 
is slightly out of focus, displaces towards it due to the 
increasing influence of other states. 




FIG. 4. Contour plot of p c (r,uj) for 8 = 50 meV (top left), 
Ap c (r, ui) for 5 = 20 meV (top right) and Ap c for 5 = 50 meV 
(bottom left) , with lo — 10 meV. For comparison, — 1 9242 (r") | 2 
is displayed at the bottom right. 

The analysis of the decoherence effects point out the 
conditions to improve the signal at the mirage, in addi- 
tion to control of the confinement. Keeping for simplicity 
Ri on the major axis and the mirage point at aRi = —R4, 
one should find a wave function ifi(r) with large ampli- 
tude at some point Ri, with the condition that all wave 
functions <fij{r) with opposite parity under a than ipi(r) 
and \ej — et\ < 0.1 eV have small amplitude at ±Ri. The 
size of the ellipse is then changed keeping e constant to 
set ei near sf and the impurity is positioned at Ri. It 



turned out to be difficult to find an ideal case for condi- 
tions close to the experimental ones. However, inspection 
of several tenths of wave functions suggests that if a and b 
are reduced by a factor 1.104 (so that a becomes 64.6 A) 
in order that state 35 reaches the Fermi energy, and the 
impurity is placed at a distance of 0.4a from the center of 
the ellipse, the feature observed at aRi should be larger 
than that of the original experiment. Keeping other pa- 
rameters as in Fig. 3 (b), the resulting Ap c (±Ri,uj) is 
shown in Fig. 5. The signal at the mirage is in fact 
nearly two times larger than in the previous case. The 
partial loss of coherence is here due to the states 28 and 
42, which are even under a, while 35 is odd. This result 
seems hard to explain in terms of billiard scattering. 
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FIG. 5. Ap c (r, ui) as a function of uo at two positions 
r = (a;,0), for an ellipse with £35 = ep, 8 = 50 meV, and 
the impurity positioned at Ri = (— 0.4a,0). 

In order to investigate if non-magnetic impurities like 
S or Si might lead to a mirage effect, I have repeated 
some calculations for U — and tuning Ed so that /3d(w) 
is peaked near ep- In agreement with Ref. Q, a mirage 
effect roughly consistent with experiment is possible only 
if the Vj are reduced by a factor ~ 5, so that the hopping 
£' ~ 0.1 eV. This seems unrealistic since replacement of 
3d by s or p orbitals is expected to increase the hopping 
and involve more sites on the surface. || In addition, the 
property of the Kondo resonance to tune itself near ep is 
lost. 

In summary, I have shown that the level width of the 
conduction states in a quantum corral, determined by 
the quality of the confinement, controls the shape of the 
observed Kondo feature and the intensity of the mirage 
effect. The nature of the decoherence effect is clarified 
and the conditions to observe an enhanced effect with 
the impurity out of the foci are stated. Replacement of 
the Kondo resonance by a phenomenological form or a 
low-energy phase shift might be valid for poor confinc- 
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ment, but for the general description, a genuine many- 
body treatment is needed. To my knowledge, this is the 
first calculation of this type for the quantum mirage. 
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side an ellipse, (h]] Helpful discussions with C.R. Proetto, 
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edged. This work benefitted from PICT 03-00121-02153 
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